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Abstract 

Superradiance of nuclear spins is considered, when the nuclei interact via 
hyperfine forces with electrons of a ferromagnet. The consideration is based on 
a microscopic model. If the sample, coupled with a resonant electric circuit, 
possesses electronic magnetization, then the electron subsystem plays the role of 
an additional effective resonator, by enhancing the coupling between nuclear spins 
and the resonant circuit. Radiation power can be increased by three orders, while 
the radiation time of a superradiance burst can become three times shorter. In 
the presence of dynamic nuclear polarization, the regime of pulsing superradiance 
can occur. 



1 Invited report at the International Workshop on Laser Physics, Bordeaux (2000). 
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1 Introduction 



The possibility of self-organized nuclear spin superradiance was predicted by Bloem- 
bergen and Pound [1] and observed in a series of experiments [2-5]. A microscopic 
theory of this phenomenon was developed in Refs. [6-8], the results being in good 
agreement with experiments as well as with computer simulation [9-12]. 

The experiments [2-5] have dealt with proton spins inside paramagnets, such as 
propanediol C 3 H 8 2 , butanol C 4 H 9 OH, and ammonia NH 3 . When nuclei are incorpo- 
rated in a ferromagnet, their spins interact with those of electrons by means of hyperfine 
forces. If the electron subsystem possesses long-range magnetic order, this can essen- 
tially influence nuclear spin dynamics [13]. The aim of this report is to describe how 
the arising electron ferromagnetism influences nuclear spin superradiance. 

The description of nonlinear spin dynamics is based on the scale separation approach 
[7,8,14]. Because of the importance of this approach for obtaining a detailed picture of 
spin evolution, we find it pertinent to briefly sketch in this introduction the main points 
of that technique as applied to spin systems. The basic parts of the scale separation 
approach [7,8,14] are: (i) Short-range stochastic quantization; (ii) Classification of 
relative quasi-invariants; and (iii) Generalized averaging technique. 

Short-range stochastic quantization makes it possible to derive a closed set of evolu- 
tion equations for a given statistical system. This is achieved by decoupling correlators 
in such a way that takes into account short-range incoherent effects. For illustration, 
let us consider spin operators Sj labelled by a site index i = 1,2, ... ,N. The evo- 
lution equations for the averages of these operators involve, as is known, binary spin 
correlators. To render the system of equations closed, one could decouple the binary 
correlators into the products of spin averages. This is what is called the mean-field 
approximation or semiclassical approximation. Such approximations take into account 
only long-range correlations but completely ignore short-range effects. The latter, how- 
ever, can be principally important for the correct description of evolution. To take into 
account short-range correlations, it is necessary to resort to a more elaborate decoupling 
of spin correlators, than the mean-field one. 

Let us present a binary spin correlator < SfS 1 ^ >, with i ^ j, in the form 



Here the factors < S" > are associated with long-range effects, in many cases permit- 
ting the usage of the uniform approximation 



The terms 5Sf describe short-range effects. For this purpose, these terms are treated 
as random variables modelling local spin fluctuations. To concretize the choice of the 
random variables, one has either to introduce a related distribution or to define the 
first moments of the random-variable products. If 8Sf are assumed to be Gaussian 
variables, then one has to define just the first two moments. Denoting the averaging 
over the random variables 5S^ as <C ... we set 



< 



S?Sf > = < S« >< Sf > + <S°> 5Sf+ < Sf > 5S? . 



<Sf> = - 




<6S? 



> = . 
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The binary correlators <C 5Sf5S^ ^> are to be defined according to the nature of the 
stochastic variables 5Sf describing local spin fluctuations [15]. For instance, we may 
set 

« 6S?6Sf » = i S(S + 1)5^ , 

where S is a spin value. 

The evolution equations for the spin components are obtained by means of the 
Heisenberg equations. Averaging these equations, with a statistical operator p(0), one 
gets 

d % 



<Sf> = - 



< 



S?, H 



> 



dt 1 h 

where H is the system Hamiltonian. The right-hand side of this equation contains 
binary spin correlators which are to be presented according to the procedure of short- 
range stochastic quantization described above. The resulting equation includes the 
stochastic terms 8Sf , hence, this is a stochastic differential equation. In this way, we 
obtain a closed system of equations although the price for this is that the evolution 
equations become stochastic. The random variables are associated with short-range, 
quantum, incoherent effects. That is why the procedure of deriving such equations can 
be called short-range stochastic quantization. 

The notion of relative quasi-invariants [15] is introduced for solving nonlinear sys- 
tems of stochastic equations in partial derivatives. Let us have a set of functions 
fn — /n( x , depending on a collection x = {x±,X2, ■ ■ ■} of variables x«, on a set 
of stochastic variables, £, and on an ensemble e = {£1,62, •••} of small parameters, 
\ek\ <C 1. The collection x can contain, e.g., spatial and a temporal variables. For 
each variable Xi, we define a variation length Li as a characteristic distance at which 
the function f n changes essentially. If the function f n is periodic with respect to x iy 
then the variation length Lj is the distance between a maximum and a minimum of 
f n for varying Xi, that is, Li is a half-period. For a nonperiodic function, the variation 
length Li can be defined as the linear size of the considered system with respect to the 
variable Xi. In particular, Lj can be infinite. Let us introduce the notation 

\/\i= I n 

in which the integration with respect to Xj is over the interval [0, Lj], and <C . . . ^> 



implies a stochastic averaging over £. 
respect to f n if, and only if, 



A function f m is an rr r quasi-invariant with 



lim 



_d_ 

dxi 



o, 



lim 

e^O 



_d_ 

dxi 



fn 



7^0, 



where the limit e — > denotes that either all Ek — ■> or that there exists a subset {Ek\ 
of the set e, such that the above limits hold true for all Ek from this subset. 

The notion of relative quasi-invariants makes it possible to classify the solutions 
of evolution equations as fast or slow with respect to each other and to generalize 
the averaging technique [16] to nonlinear systems of stochastic equations in partial 
derivatives. 
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2 Electron-Nuclear Hamiltonian 



We consider a Hamiltonian of the general form 

H = H e + H n + H en , (1) 
describing interacting electrons and nuclei. In the electron Hamiltonian 

H e = — - E "A? S? • Sj — // e ^2 B ■ Sj , (2) 

i^j i 

where is an exchange interaction; Sj, an electron spin operator; \i e = g e [J>B, with g e 
being the electronic gyromagnetic ratio and /i B , the Bohr magneton; B is a magnetic 
field. The nuclear Hamiltonian 

H n = \ll E Cglflf-Hn E B I, (3) 

i^j a/3 i 



contains the dipole interactions 



Cf = & (5 a , - 3< ng) , 
ij 



where n n = g n ^N, with g n being the nuclear gyromagnetic ratio; fijy is the nuclear 
magneton; = | r^- 1 , = r j — r,-, n^- = r^/r^; and Ij is a nuclear spin operator. The 
electron-nuclear interactions are described by a hyperfine Hamiltonian 

^^ESrU^EE^^^- ( 4 ) 

i ij^j a/3 

in which A is an isotropic contact interaction and 

A a/3 _ l^ePn fr % a ri /3\ 

' ij 

is a dipole interaction between electron and nuclear spins. The total field 

B = H e z + H x e x , Hi = H a + H (5) 

consists of an external field H Q and the field Hi, in which H a is an effective field of a 
transverse magnetocrystalline anisotropy and H is a resonator feedback field. 

The sample is coupled to a resonant electric circuit of inductance L, capacity C, 
and resistance R defining the natural frequency u, ringing time 73, and quality factor 
Q by the relations 

1 UJ ujL 

U1 = 7W' 73 = ^' Q = 1T- 

The electric current in the resonator coil is induced by the motion of the transverse 
magnetization 

M * = y E (»e < S! > +fi n < Jf >) , (6) 

i 



4 



where V is the sample volume. The resonator feedback field is described by the Kirch- 
hoff equation 

— + 2 l3 H + uo 2 / H{r) dr + Aur] = , (7) 
at Jo or 

where 77 is a coil filling factor. 

To write down the evolution equations, it is convenient to pass to the ladder spin 
operators 

S j = S j ± lS i > I i = T l ± lI j ■ 
Then the electron spin Hamiltonian (2) takes the form 

He = -\Y. Ja(st Sj + sts^-^Y.\H,st + l -H 1 (si + sr) 

The nuclear Hamiltonian (3) becomes 

Hn = \Y. (ey It I' + Oy If // + 6y /+ if + Jf J," + 2Cy /+ //+ 



where 



n . . = r zz p . = — ( r xx 4- r yy \ 

a i3 — U ij ' e *J — 2 V »•?' *■? / ' 

The hyperfine Hamiltonian (4) transforms to 

4„ = ^ x: ir +s-it+ 25? /?; 



with 



+ 2 1 ^(S+Ir + Sri;) + ^ S? I* + ^ S+ 1+ + S~ I~+ 

(S? 1; + s? 1+) + (4)* (.s, // • .s; / ; ) 



(8) 



(9) 



(10) 
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3 Evolution Equations 

From the Heisenberg equations for the spin operators we derive the equations of motion 
for the ladder and longitudinal spins. For the ladder electron spins, we have 



dt 



= -VL e (H Sr-H 1 S?)- y E Jy (Sr Sj - S? Sr) + 



+A(sr ii-sii-) + \ £ sr i; - 4, s; i--2V ij s?i+ + 

3\+i) 

+ Cj • S i If + (c- ■)* S i L — 2c- • S? I z 



where and in what follows we set, for the simplicity of notations, the Plank constant h = 
1. An equation for St is obtained from Eq. (11) by means of Hermitian conjugation. 
For the longitudinal electron spins, we get 

+\ E Jv(Sr S} -St Sj) + \A{st ir -Sr + 

+\ E \\ 4 (st ij - s- 1+) + ^ st it - 



- %)* sr i- + c' i3 st /; - (4)* s r i 

For the ladder nuclear spins, we find 

. di; 



(12) 



dt 



-ix n (H I' - H x It) + 



+ E («« K Ij ~ I! IJ -26y/?/; + Cy I' It + c?./r/r - 2cy /?//) + 



3&i) 



+A(lr S? - I? Sr) + \ £ Ir S* - ^ I* S~ - 2b> I* 5+ + 



2 • 
J i(^) 



+ ^/r57 + (^)*/r57-2c^s- 

Finally, for the longitudinal nuclear spins, we obtain 



(13) 



dJl 
dt 



^ nHl (l--lt) + 



+ V [- e- • (7+ Jr - I~ + 5- • J + /+ - 6*. J.~ J~ + c- ■ J + J 2 - c*- J-~ J 2 



3{+i) 



+ 



+lA(its--i-st) + 1 - E \\<3 {itsj-Ksfj^itsf- 
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- W K S J + c 'ij !f S j - W K S]\ . (14) 

In the calculations below, we shall employ the following properties of the dipole 
interactions: 

E A? = £ c- = o, E <= E c<f = 0. 

From the first of these properties, it follows that 

I I ' 

The second of the summation properties above is, strictly speaking, approximate being 
valid up to boundary effects. 

The Zeeman frequencies for electrons and nuclei are 

U e = fl e H , U n = ll n H , (15) 

respectively. The characteristic wavelengths 2n/u e and 2n/u n are much larger than the 
mean distance between spins. Therefore, for the statistical averages of spin operators 
of electrons, 

x = < S~ > , z = < S? > , (16) 

and nuclei, 

u=<I~>, s = <I->, (17) 

we may use the uniform approximation. At the same time, local spin fluctuations, 
disturbing space uniformity, will be taken into account by means of the short-range 
stochastic quantization explained in the Introduction. Realizing this procedure, we 
come to the following expressions of local random fields: 



£o = \ E fan fa? - 6Sf) + a'^I* + 4,61+ + (c^)* SI- 

^ = 1 E [2^ {SST ~ SS~) + e ' tJ 5IJ + 2b> J 51+ + 2^ 51* 
^ = \H [3% 51* + 2 Cij 51+ + 24 SI' + ^ 5S* + 4 5S+ + (^)* 5S j 

V = \ E [ 4 % SI+ + 4cy 5/; + 4 55/ + 26(,. 5S+ + 24 55?] . (18) 

Introduce also the anisotropy frequencies 

a e = fi e H a , a n = /i„ # a . (19) 

Averaging the Heisenberg equations (11) to (14), we employ the short-range stochastic 
quantization, include into the equations the transverse and longitudinal spin relaxation 
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parameters, and envolve the notations (15) to (19). This results in the evolution 
equations for the electron transverse spin average, 

dx 

— = i{u: e - is-( + 272)^ - i(a e - A u - £ + fx e H)z , (20) 

electron longitudinal spin average, 

dz 1 ' 

— = - (a e - A u - i + fi e H) x* - -{a e - Au* - f + p e H)x - lx {z - a) , (21) 

at 2 2 

where a is a stationary electron single-spin polarization, for the electron transverse 
modulus of spin squared, 



d\x 



2 

- = -2^ 2 \x\ 2 + i(a e - A u* - C + fi e H)zx - i(a e -Au-£ + fjL e H)zx* , (22) 



dt 

for the transverse nuclear spin average, 
du 

— = i(u n - A z - (p + iT 2 )u - i(a n - A x - if + \i n H)s , (23) 
dt 

longitudinal nuclear spin average, 
ds 1 % 

— = - (a n - A x - up + fi n H) u* - - (a n - A x* - cp* + fi n H) u - Y^s - Q , (24) 

where ( is a stationary nuclear spin polarization, and, finally, for the nuclear transverse 
modulus of spin squared, 

d\u\ 



dt 



= — 2r 2 \u\ + i (a n — A x* — p* + Hn H) su — i (a n — A x — p + //„ H) su* . 

(25) 

The evolution equations (20) to (25) form a set of nonlinear stochastic differential 
equations, which are also to be complemented by the Kirchhoff equation (7), where the 
magnetization (6) has to be understood as the total average 

M x = X - < [fi e p e [x + x*) + n n Pn {u + u*)\ > , (26) 

with p e and p n being the electron and nuclear densities, respectively. It is convenient 
to rewrite the feedback equation (7) in the form 

H = -Ati V f G{t - r) dM x (r) , (27) 
Jo 

with the Green function 

G(t) = (cosuj 3 t — — sincj 3 t ) e~ 73 * , uj 3 = J uo 2 — 7! . 
V uj 3 J 

To make the set of the evolution equations completely defined, we need to concretize 
the random variables (18). For this purpose, we treat these variables as Gaussian, with 
the stochastic averages 

<£ > = <e> = 0, <£ £> = 0, <£ 2 > = < |£| 2 > = 7 2 (28) 
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for the electron-spin local fluctuations, and with 

<y?o> = <V 5 > = 0, <v ? oV 5 > = 0, <y?o> = < W? > = ?l (29) 

for the nuclear-spin local fluctuations. The quantities 7* and are inhomogeneous 
widths, for which, according to the definition (18), one may write 

2 -p2 1 -p2 -p2 -p2 1 -p2 f^fl^ 

7* 1 ee ~r 1 en i 1 * L en * L nn ) \ ou / 

where the corresponding terms are due to electron-electron, electron-nuclear, and to 
nuclei local interactions. 



4 Small Parameters 

To simplify the evolution equations (20) to (25), let us take into account the existence of 
several small parameters. First, the longitudinal and transverse relaxation parameters 
are assumed to be small as compared to the Zeeman frequencies, so that 

— < 1 , — < 1 , — < 1 , — < 1 . (31) 

UJ e U e UJ n UJ n 

Inhomogeneous broadening is supposed also to be week, 

^ « 1 , — « 1 . (32) 
The same concerns the magnetocrystalline anisotropy, 

- « 1 , — « 1 . (33) 

The energy of the spin interaction with the resonator feedback field is much weeker 
than the corresponding Zeeman frequencies, which implies that 

2 

Pe /^e ^ Pra Me pn ^ ^ 

PePePn < ± ^ PnJ^ < j ^ 

This is equivalent to saying that \ji e H\ <C u e and \/i n H\ <C uj n . The contact hyperfine 
interaction is smaller than the electron Zeeman frequency, 

— < 1 , (35) 

but A can be comparable or even much larger than uj n . The resonator is of good quality, 
i.e. its quality factor is high, Q S> 1, which, because of the relation 73 = u/2Q, means 
that the resonator ringing time is small, 

— < 1 . (36) 

UJ 



Since the nuclear magneton is three orders smaller than the Bohr magneton, one has 



And also, one usually has 



^ « 1 . (37) 



— < 1 , — < 1 . (38) 
7i 72 



Due to the inequality (37), the Zeeman frequencies are related as 



— < 1 • (39) 



Therefore, the resonant electric circuit can be tuned either to u n or u e . 

In general, the dynamics of electron spins is similar to that of nuclear spins. The 
main difference is that nuclear spins, owing to the inequality (37), or to that (39), 
weekly influence the behaviour of electron spins. While, to the contrary, electronic 
spins, can essentially influence spin dynamics. For instance, the effective electronic 
spin-resonance frequency, as is seen from Eq. (20), is shifted as u e — As, as a result of 
the hyperfme interaction. But due to the inequality (35), this shift is very small, and 
can be neglected. To the contrary, the effective nuclear magnetic- resonance frequency, 
as follows from Eq. (23), is shifted as u n — Az, which is a kind of the dynamical 
frequency shift [17,18]. In the presence of the long-range magnetic order in the electron 
subsystem, the nuclear magnetic-resonance frequency becomes 

u N = co n - A a , (40) 

where a is a stationary electron magnetization. Since electrons, especially those pos- 
sessing magnetic order, can essentially influence the evolution of nuclear spins, but not 
conversely, the nuclear spin dynamics exhibits more varieties and is more interesting 
than that of electronic spins. Therefore, we shall concentrate on the nuclear spin dy- 
namics, and will imply in what follows that the resonator is tuned to the magnetic 
resonance frequency (40), so that the quasiresonance condition 



U! N 



<1, A N = u-u N (41) 



holds true. 

Experimentally, the radiation intensity of moving spins can be observed through 
measuring the current power 

9 9 V H 2 

P = RJ 2 , J — 



Att i] L 



where R and L are resistivity and inductance of the resonant electric circuit. Thus, we 
have 



7s V u2 



27T T) 



H 2 . (42) 
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One may analyze how the arising magnetic order changes the power (42) or its time 
average P, with an averaging over fast oscillations. Considering P = P(cr) as a function 
of electron magnetization, one can study the relative difference 

In order to estimate the values (42) and (43), one has to calculate the resonator field 
(27), which, under the inequality (36), envolves the Green function 

G(t) = cos(uj t) e~ 73 * . 

Estimates show [13] that the relative quantity (43) is 

SP( a ) ~ P^eAa (44) 

Pn Pn 

Hence, if p e ~ p n and un ~ A, the power (42) can be increased by three orders when 
a ~ 1, since then 

SP( a ) ~ ^ a > 1 , (45) 

Pn 

which follows from inequality (37). This is the enhancement effect caused by the 
electron magnetization. 

5 Time Evolution 

The existence of small parameters, discussed in the previous section, allows us to clas- 
sify the solutions to the evolution equations (20) to (35) onto fast or slow as compared 
to each other. Thus, among all functions, x is the fastest one. In other words, all solu- 
tions z, \x\ 2 , u, s, and \u\ 2 are temporal quasi-invariants with respect to x. Dynamics 
of nuclear spins is described by equations (23) to (25). Averaging these equations over 
the period T e = 2^/uo e of the fastest oscillations, related to x, and setting 



1 rt+T e I rt+T e 

/ x(t) dr = , — / z(t) dr = a , 



we have 



du ., . ., — 

— = i{uj N - ipo + iT 2 ) u - i{a n -<p + ti n H)s, 
at 

^ = l -{a n - ip + n n H) u* - ^(a n - ip* + /j, n H) u - T^s - () , 
and similarly for \u\ 2 , where 

T e It 

From here, the functions s and \u\ 2 are quasi-invariants with respect to u. 



rt+T e 
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Following further the averaging technique, one can derive the equations for the 
slowest functions s and \u\ 2 . The latter equations allow analytical solutions in the case 
of short time f C T 2 * = T^ 1 . This analysis can be found in Refs. [6-8,13]. The appear- 
ance of pure spin superradiance depends on the value of an affective coupling between 
the nuclear spins and the resonator and also on the level of the initial spin polarization. 
The latter can be varied in a large diapason ranging from zero to practically 100% [19]. 

For larger times, when t ~ T 2 * or t ^> T 2 *, it is necessary to accurately consider the 
inhomogeneous broadening. Then the analytical solutions of the evolution equations 
is not available and one has to resort to numerical calculations. When the polarization 
of nuclear spins is supported by the procedure of dynamical nuclear polarization, so 
that the parameter ( < in Eq. (24), then the regime of pulsing spin superradiance 
can arise [20]. 

To illustrate the general behaviour of nuclear spin dynamics, we present the results 
of numerical calculations. Figure 1 shows the regime of weak superradiance [8] without 
dynamical nuclear polarization, while Figs. 2 to 8 demonstrate the regime of pulsing 
spin superradiance occurring in the presence of dynamical polarization. Time is every- 
where measured in units of T^ 1 . In Figs. 2 to 8, the effective coupling parameter is 
set g = 10, and the pumping parameter is ( = —0.5. The function w(t) = \u(t)\ 2 is 
proportional to the current power. This is why we concentrate our attention on this 
behaviour. We study the behaviour of w(t) for the varying parameter 7 = ri/r 2 and 
for different initial conditions w = w(0) and s = s(0). To simplify calculations, the 
resonator feedback field was treated in an approximation that neglects inhomogeneous 
broadening, because of which the presented figures should be considered as only a 
qualitative illustration. But we do hope that the regime of pulsing superradiance will 
survive in a more elaborate treatment that is in progress. 
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Figure Captions 

Fig. 1. The dipole radiation intensity I(t) of nuclear spins, the coherence coefficient 
Ccohit), and the polarization p z (t) = —s(t) as functions of time. The radiation intensity 
is given in arbitrary units, and time is measured in units of T^ 1 - This figure shows the 
regime of weak superradiance [8]. 

Fig. 2. The regime of pulsing spin superradiance for 7 = 0.01, Wq = 10~ 6 , and 
s = -0.1: (a) w(t); (b) s(t). 

Fig. 3. The function w(t) versus time for 7 = 0.001, w = 10~ 6 , and s = —0.1. 

Fig. 4. The solution w(t) as a function of time under the initial conditions Wq = 0.1 
and s — —0.25 for varying 7: 7 = 1 (solid line) and 7 = 0.5 (dashed line). 

Fig. 5. The function w(t) at wo = 0.5 and s = 0.5 for 7 = 1 (solid line) and 
7 = 0.5 (dashed line). 

Fig. 6. Pulsing spin superadiance for 7 = 0.01, w = 0.01, and s = —0.1. 

Fig. 7. Pulsing superadiance in the case of the parameters 7 = 0.01, w = 0.001, 
and s — —0.5. 

Fig. 8. The regime of pulsing spin superradiance for 7 = 0.01, Wo = 0.1, and 
s = -0.1. 
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